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Abstract 

Experimental data on the number of neutrons born in the heavy water targets 
of the large neutrino detectors are used to set the limit on the proton life-time 
| independently on decay mode through the reaction d — > n+1. The best up-to-date 

limit t p > 4xl0 23 yr with 95% C.L. is derived from the measurements with D2O 
target (mass 267 kg) installed near the Bugey reactor. This value can be improved 
by six orders of magnitude with future data accumulated with the SNO detector 
containing 1000 t of DoO. 
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1 Introduction 

o : 

The baryon (B) and lepton (L) numbers are absolutely conserved in the Standard Model 
D ■ (SM). However, many extensions of the SM, in particular, grand unified theories incorpo- 
rate B and L violating interactions, since in the modern gauge theories conservation of 
baryon (lepton) charge is considered as approximate law due to absence of any underlying 
symmetry principle behind it, unlike the gauge invariance in electrodynamics which guar- 
antees the massless of photon and absolute conservation of the electric charge. Therefore, 
^ ■ it is quite natural to suppose the decay of protons and neutrons bounded in nuclei. The 
processes with AB = 1, AB = 2, A(B — L) = 0, A(B — L) = 2 have been discussed in 
literature (see e.g. |], 0, § and references therein), while the disappearance of nucleons 
(or decay into "nothing") has been addressed in connection with possible existence of 
extra dimensions @, §, |J . 

Stimulated by theoretical predictions, nucleon instability has been searched for in 
many underground experiments with the help of massive detectors such as 1MB, Frejus, 
Kamiokande, SuperKamiokande and others (for experimental activity see || [7], ^] and 
refs. therein). About 90 decay modes have been investigated; however, no evidence for 
the nucleons decay has been found. A complete summary of the experimental results is 
given in the Review of Particle Physics |J. For the modes in which the nucleon decays 
to particles strongly or electromagneticaliy interacting in the detector's sensitive volume, 
the obtained life-time limits are in the range of 10 30 — 10 33 yr, while for decays to only 
weakly interacting products (neutrinos) the bounds are up to 10 orders of magnitude 
lower P, However, because it is not known a priori which mode of proton decay 



(from 90 ones listed in ||) is preferable, the limits on the proton decay independent on 
the channel are very important. 
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Three approaches were used to establish such limits: 

(1) In ref. ]TTJ the bound r(p — >?) > 1.3xl0 23 yr was determined^ on the basis of 
the limit for the branching ratio of 232 Th spontaneous fission /3gp- It was assumed that 
the parent 232 Th nucleus is destroyed by the strongly or electromagnetically interacting 
particles emitted in the proton decay or, in case of proton's disappearance (or p decay 
into neutrinos) by the subsequent nuclear deexcitation process. Using the present-day 
data on the 232 Th half-life T 1/2 = 1.405xl0 10 yr and limit (3 S f < 1.8xl(r 9 %, we can 



recalculate the value of |Tl] as r(p — >?) > l.OxlO 23 yr. 

(2) In ref. |13| the limit r(p — >?) > 3xl0 23 yr was established by searching for 
neutrons born in liquid scintillator, enriched in deuterium, as result of proton decay in 
deuterium (d — > n+?). 

(3) In ref. |14] the limit r(p — > 3u) > 7.4xl0 24 yr was determinedP] on the basis of 
geochemical measurements with Te ore by looking for the possible daughter nuclides ( 130 Te 



129- 



Xe), while in refs. Hl5| , |16|| the bound r(p — ■> 3v) > 1.1x10 yr was achieved 



by the radiochemical measurements with 1710 kg of potassium acetate KC2H3O2 placed 
deep underground ( 39 K — > 37 Ar). In the experiments (2) and (3) both the baryon 

number and the electric charge would be not conserved; nevertheless authors suggested 
that "experimenter would be wise not to exclude such processes from consideration a 

priori" 
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The limits |L4|, [T^, [L6| usually are quoted as "independent on channel" ||, however 
it is evident that they are valid only for the proton decay into invisible channels (or 
disappearance), in which the parent nucleus is not fully destroyed (like 232 Th in the 
experiment [llj). At the same time, bound on the proton decay from the deuterium 
disintegration requires the less stringent hypothesis on the stability of daughter nuclear 
system and, hence, it is less model dependent. Such a limit was established in 1970 
13 1 and is equal r(p — ►?) > 3xl0 23 yr at 68% CL.Q This value can be improved by 



using the data from the modern neutrino experiments with heavy water well shielded 
against cosmic rays and natural radioactivity. With this aim, in present paper we analyze 
the measurements of ref. |T7| with the 267 kg D2O target installed at Reactor 5 of the 
Centrale Nucleaire de Bugey (France). Further, we show that obtained limit r(p — >?) can 
be highly improved with the SNO (Sudbury Neutrino Observatory) large volume detector 
|18| developed mainly for the Solar neutrino investigations and containing 1000 t of D2O. 



2 D 2 experiment at the Bugey reactor 



The experiment [|17]] was aimed to measure the cross sections for the disintegration of 
deuteron by low-energy electron antineutrinos from nuclear reactor through reactions 
V e + d—>-T7 e + n + p (neutral currents) and V e + d~^e + +n + n (charged currents). 
Events were recognized by the neutrons they produced. The detector was located on 
the depth of 25 meters of water equivalent (mwe) at 18.5 m distance from the center of 



1 We recalculated the value quoted in |TT| t(N — »?) > 3 x 10 23 yr (given for 232 particles: 142 neutrons 
and 90 protons) for 90 protons which should be taken here into consideration (JV is p or n). 

2 The value t(N — > 3^) > 1.6xl0 25 yr quoted in juj as given for 52 particles (28 neutrons and 24 
protons) was recalculated for 24 protons. 

3 Because ref. jl3| is not the source of easy access, and in ref. ||, where this limit is quoted, there is 
no indication for confidence level, we suppose that it was established with 68% C.L. 
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the Reactor 5 core at the Bugey site. The cylindrical target tank, containing 267 kg of 
99.85% pure D 2 0, was surrounded by layers of lead (10 cm) and cadmium (1 mm) to 
absorb thermal neutrons from external surroundings. The tank with D 2 and Pb-Cd 
shield was inserted in the center of a large liquid scintillator detector (based on mineral 
oil) which served as (inner) cosmic ray veto detector. Subsequent layer of lead 10 cm thick 
was aimed to reduce the flux of external 7 quanta with energies E 1 > 2.23 MeV which can 
photodisintegrate the deuterons creating the background events. However this shielding 
itself was a significant source of neutrons in the target detector: they were created due to 
interaction of cosmic rays with Pb. To suppress this background, an additional layer of 
cosmic ray veto detectors was installed outside the Pb shielding. The outer veto reduced 
the neutron background in the target by a factor of near 6. Neutrons were detected by 
3 He proportional counters installed in the tank with D 2 through reaction 3 He + n — > 
3 H + p + 764 keV. Further details on the experiment can be found in [17]. 



The decay or disappearance of proton bounded in deuterium nucleus, which consists 
only of proton and neutron, will result in the appearance of free neutron: d — > n+7. Thus 
the proton life-time limit can be estimated on the basis of the neutron rate detected in 
the D 2 volume when the reactor is switched off. To calculate the limr(p — ►?) value, we 
use the formula 

lim r(p -*?) = e x N d x t/limS, (1) 

where e is the efficiency for the neutron's detection, Nd is number of deuterons, t is 
the time of measurement, and \imS is the number of proton decays which can be ex- 
cluded with a given confidence level on the basis of the neutron background (one-neutron 
events) measured in the experiment. Mean efficiency for single neutrons born isotropically 



throughout the D 2 volume was determined as e = 0.29±0.01 [pLY|] . In 267 kg of D 2 
there is Nd = 1.605 xlO 28 deuterons. Raw one-neutron rate with the reactor down is equal 
25.28±0.68 counts per day (cpd) and, corrected for software efficiency (0.444±0.005), this 
rate is 57.00 ±1.53 cpd. For very rough estimate of the p life-time (as the first approxi- 
mation) we can attribute all neutron events to proton decays and obtain the lim S value 
as 59.5 cpd at 95% CL. Then, substituting this value in the formula (1) we get limit 
r(p — ►?) > 2.1xl0 23 yr with 95% C.L.R However, it is evident that r limit derived in 
this way is very conservative because the dominant part of observed neutron rate has 
other origins rather than proton decay | ]2"0"| , |2~T| , 22" . On the other hand, if we suppose that 



all measured neutron events are belonging to background, then the excluded number of 
neutrons due to possible proton decay will be restricted only by statistical uncertainties 
in the measured neutron background. Hence, to estimate value of lim S we can use so 
called "one (two, three) a approach", in which the excluded number of effect's events is 
determined simply as square root of the number of background counts multiplied by one 
(two or three) according to the confidence level chosen (68%, 95% or 99%). This method 
gives us the sensitivity limit of the considered experiment to the proton decay. Applying 
it we get limS = 3 cpd (at 95% C.L.), which leads to the bound r(p — >?) > 4xl0 24 yr. 

Therefore, we can argue that, with the probability close to 100%, estimate of r limit is 
within interval 2xl0 23 yr - 4xl0 24 yr. In order to fix life-time limit or at least to narrow 
this interval, it is necessary to determine the contributions of different sources to the total 



4 The similar limit r(p — >?) > 1.9xl0 23 yr with 95% C.L. can be derived from other neutrino deuteron 
experiment at Krasnoyarsk (Russia) nuclear reactor ]19[. 
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Table 1: Main characteristics of the neutrino experiments at nuclear reactor 



Experiment 


Bugey-1995 

El 


Bugey-1999 



Palo Verde 

n 


Chooz 

n 


Depth 


40 mwe 


25 mwe 


32 mwe 


300 mwe 


Distance 


15 m 


18.5 m 


^0.9 km 


wl km 


Detector type 


Li-loaded 
scintillator 


D 2 target + 
3 He counters 


Gd-loaded 
scintillator 


Gd-loaded 
scintillator 


Detector mass 


^0.6 t 


267 kg 


11 t 


5 t 


Neutron 
background 


^100 cpd/t 


^100 cpd/t 


2.2 cpd/t 


0.24 cpd/t 



neutron rate observed. As it was already mentioned, the nature and origins of neutron 
background in neutrino experiments at nuclear reactors are well known and understood 
(see, for example, refs. 



20, 21, 22 



I). The main sources are: (i) interaction of cosmic 
muons (escaped an active veto system) with the detector, passive shield and surrounding 
materials; (ii) photodisintegration of the deuteron by 7 quanta (with E 1 > 2.23 MeV), 
originated from the radioactive contamination of the detector materials and shield, as well 
as from environment pollution; (iii) residual (and non-eliminated by the shield) neutron 
background at the nuclear reactor site. 

Before coming to details, we would like to remind that crucial characteristics of any 
neutrino experiment at reactor are the depth of its location and distance from the reactor 
core [pO , pl| , |22 , p3| . Let us prove this statement by Table 1 with parameters of the most 
advanced experiments and by two short citations, from ref. f20| : "A striking difference 
between the two experiments is the amount of overburden, which may be viewed as the 
main factor responsible for how the experiments compare on detector design, event rate, 
and signal-to-background", and from ref. [F2TJ: "... a detector should be located sufficiently 
deep underground to reduce the flux of cosmic muons - the main source of background in 
experiments of this type". 

It is clear from the Table 1 that neutron background of different experiments is de- 
creased as the depth of their location and distance from the reactor are enlarged. For 
example, in the Chooz experiment with liquid scintillator the background was reduced 
roughly by factor 500 as compared with that of Bugey ones fl7|, |24| because the Chooz set 
up was placed deep underground (300 mwe) and ~1 km away from nuclear reactor ||26|| . 
We recall that Bugey set ups |17|, |4j] were located only 25-40 mwe overburden (which 
allows to remove the hadronic component of cosmic rays but is not enough to reduce 
the muon flux significantly) and at 15-18 m distance from the reactor core. Thus, the 
dominant part of neutron background in [[H], is associated with the reactor site and 
muon flux. Indeed, as it was proved by the detail simulation and careful analysis of neu- 
tron background in reactor-off periods of the experiment [2~1|, the 67±3% of neutron rate 
(measured at 15 m distance) are attributed to known origins (see Table 6 in ref. |[24|| ). 



On the basis of comparison of different experiments presented in Table 1, and taking into 
account results of background analysis pl |, we can make semi-quantitative and conser- 
vative estimation that at least 50% of one-neutron events measured in [17] are caused by 



4 



the mentioned sources (i) - (iii). Hence, attributing remaining part of one-neutron rate 
to other unknown background origins, we can accept its value as the excluded number of 
proton's decays (limS 1 = 30 cpd). Finally, substituting this number in the formula (1) we 
obtain 

r(p ->?) > 4xl0 23 yr with 95% C.L., 



which is higher than previous limit 13 . 



3 Expected improvements with the SNO solar neu- 
trino detector 

The Sudbury Neutrino Observatory (SNO) |TE| is an unique large Cherenkov detector 
constructed with an emphasis on the study of Solar neutrinos. The detector, containing 
1000 t of 99.917% isotopically pure heavy water, is located in the INCO Creighton nickel 
mine near Sudbury, Ontario, on the depth of 2039 m (near 6000 mwe); this reduces the 
muon flux to 70 muons per day in the detector area. Particular attention is payed to 
minimization of radioactive backgrounds. Near 7000 t of ultra-pure light water shield the 
central D 2 detector from natural radioactivity from the laboratory walls. All components 
of the detector are made of selected materials with low radioactivity contamination. 

Solar neutrinos will be detected through the following reactions with electrons and 
deuterons: Vi + e~ — > Ui + e~ (elastic scattering; i = e,/i,r), v e + d —> e~ +p + p (charged 
current absorption) and z/j + rf— > h>i + n + p (neutral current disintegration of deuteron). 
Near 9600 photomultiplier tubes are used to observe the Cherenkov light produced on 
the D 2 volume by high energy products. Neutrons released in d disintegration will be 
detected by neutron capture on deuterons in pure D 2 0, or by capture on 35 C1 by dissolving 
MgCl salt in the heavy water, or by capture on 3 He using proportional counters. Further 
details can be found in Jl8| . 



Extensive Monte Carlo simulations were performed to predict response functions and 
numbers of expected events due to interaction of the detector with Solar neutrinos, natural 
radioactivity of various detector components, cosmogenic activities, capture of neutrons, 



(a, fry), (a, wy) reactions outside the SNO detector, etc. (f.e., see p7fl). Expected number 
of neutrons from all sources in the D 2 volume is calculated as ^5x 10 3 during 1 yr period 
of exposition, with main contribution from the Solar neutrinos. Efficiency for n detection 
is equal 83% for n capture on 35 C1 |1| |27 |. 



Using these unique features of the SNO detector (super-low background, large amount 
of D 2 and high sensitivity to neutrons), the limit on the proton decay independent on 
channel can be highly improved. Again for very rough estimate of the p life-time we can 
conservatively attribute all neutrons in the D 2 volume to proton decays and accept it 
as the excluded value limS* = 5xl0 3 counts. Then, substituting in the formula (1) values 
of efficiency e = 0.83, measuring time t — 1 yr, number of deuterons iV d = 6 x 10 31 and 
limS* = 5xl0 3 counts, we receive 

r{p ->?) > lxlO 28 yr, 
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which is about five orders of magnitude higher than present-day limit. However this 
value can be improved further by accounting the neutron events originating from Solar 
neutrinos and high energy 7 quanta. Number of neutrons born in the D 2 volume due to 
disintegration i>i + d^i>i + n + p can be estimated independently using the information 
on the number of Solar neutrino interaction with the detector volume through neutrino- 
electron elastic scattering V{ + — > v\ + e~ . Neutrons created by high energy 7 quanta 
from natural radioactivity in the detector components can be also calculated if the levels 
of pollution of all materials are measured firmly. In that case the excluded number of 
neutrons due to possible proton decay will be restricted only by statistical uncertainties 
of the measured neutron background, i.e. we can apply "two a approach" again. It gives 
limS = 2^5000 with 95% C.L., thus the corresponding bound on the proton life-time 
would be equal to 

r(p ->?) > 4xl0 29 yr with 95% C.L., 

which can be considered as the maximal sensitivity of the SNO detector for the proton 
decay independent on channel. 

4 Conclusion 



The data of the Bugey experiment |T7|], aimed to measure the cross sections for the 
deuteron disintegration by antineutrinos from nuclear reactor, were analyzed to set the 
proton life-time limit. The obtained value r(p — >?) > 4x 10 23 yr at 95% CL. is higher than 
the limit established in the previous experiment fl3| . With the future data from the SNO 
Solar neutrino detector (containing 1000 t of heavy water) life-time limit will be improved 
up to the value r(p — >?) > 4xl0 29 yr, which is, in fact, close to the bounds established 
for the particular modes of the nucleon decays to charged or strongly interacting particles 
and would be of a great importance for many extensions of the modern gauge theories. 
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